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We evaluate the dominant superconducting pairing symmetry in some 3-dimensional cubic 
lattice structures within the third-order perturbation theory with respect to the on-site coulomb 
repulsion. We clarify whether the vertex correction term, which has a critical contribution to 
the p-wave state in the 2-dimensional systems, is also important in the 3-dimensional system. 
We finally investigate the interplay between dimensionality and the vertex correction to grasp 
common features independent of the lattice structures in the mechanism of the superconduc- 
tivity. 
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§1. Introduction 

In the strongly correlated fermion systems, the origin of superconductivity has been investigated 
with the special interest. In quasi-2-dimensional systems, such as the high-Tc cupretes^~^\ the 
d-wave superconductivity is considered to originate from the antiferromagnetic spin fluctuation. 
Such a spin fluctuation is a candidate for the mechanism of the superconductivity in the strongly 
correlated systems. Recently, on the basis of the spin fluctuation mechanism, the possibility of the 
p— and d— wave pairing states has been investigated in a variety of materials^"^^ . For example, 
Arita et al%l studied the pairing instabilities mediated by the spin fluctuation with the fluctuation 
exchange (FLEX) approximation, in the 2- and 3-dimensional Hubbard model. In their results, 
d-wave pairing is easily induced by the antiferromagnetic spin fluctuation in the high density. 
The p-wave is dominant due to the ferromagnetic spin fluctuation, in very low density. As their 
conclusion, the p-wave instability mediated by the ferromagnetic spin fluctuations is weaker than the 
d-wave instability mediated by the antiferromagnetic spin fluctuations in both 2- and 3-dimensional 
systems. Actually, the simple spin fluctuation theory does not lead to the p-wave triplet state in 
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Sr2Ru04.t-|) This indicates that the simple spin fluctuation mechanism is not the unified theory, 
but a limited one in the strongly correlated superconductors. On the other hand, the third-order 
perturbation theory (TOPT) with respect to the on-site Coulomb repulsion U can most naturally 
explain the triplet state in Sr2RuO4.0^ In this case, the p— wave triplet state is promoted by 
the wave-number dependence of the third-order vertex correction rather than the spin fluctuation 
process. 

Motivated by this fact, in the previous paperjll' we have investigated in detail the effect of the 
vertex correction in the 2-dimensional square lattice. Then, we have determined the pairing sym- 
metry in the phase diagram for electron density and next-nearest-neighborer hopping parameter. 
We have shown that the triplet state induced by the vertex correction term is dominant in the wide 
range of electron density. Thus, we have indicated that the vertex correction term has a critical 
contribution to the triplet state in the superconductivity induced by the on-site Coulomb repulsion 
in the 2-dimensional systems. This conclusion continuously presents a question whether the vertex 
correction is also important in the 3-dimensional systems. We in this paper evaluate the dominant 
pairing symmetry in some cubic lattice structures, and investigate the interplay between dimen- 
sionality and the vertex correction to grasp common features independent of lattice structures in 
the mechanism of the superconductivity. 

Before starting our study, let us discuss possibility of the magnetic instability. The strict ground 
state, actually, is determined by comparing total energies in the superconducting and the magnetic 
states. However, since our aim is not to obtain the complete phase diagram, we hereafter ignore 
the possibility of the magnetic phase transition. 

This article is organized as follows. In §2, we explain the formalism of the calculation. In §3, 
we describe the results of the calculation. We show the phase diagram and discuss the role of the 
vertex correction in detail. In §4, we summarize the interplay between dimensionality and the vertex 
correction to grasp common features independent of the lattice structures in the mechanism of the 
superconductivity, comparing with that of the simple spin fluctuation^~^\ Finally, we mention 
relation of our results to the actual materials. 

§2. Formulation 

The on-site repulsive Hubbard Hamiltonian in cubic lattice structures is given by 

'H = -ti + t2 4,aCk,a + UJ2 ni,^ni,i, (2.1) 

<i,j>,<j <i,k>,(T i 

where a is the spin index, < i^j > indicates taking summation over the nearest-neighbor sites and 
< i,k > over the next-nearest-neighbor sites. We obtain energy dispersions from the non-interacting 
part in eq. (2.1); 

Ef*^ = — 2ti(cosA;^ + cosky + coskz) + 4:t2CoskxCoskyCoskz, 
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= —8ti{coskxCoskyCoskz) + 2t2(cos(2fea;) + cos{2ky) + 003(2^2)), 

E^'-'"-^ = —4:ti{coskxCosky + coskycoskz + coskzcoskx) + 2t2(cos(2A;j;) + cos{2ky) + cos(2A;2)), (2.2) 

for the simple cubic (SC), body-centered cubic (BCC) and face-centered cubic (FCC), respectively. 
To compare the result of the cubic lattice with our previous result of the 2-dimensional square 
lattice structure (SQ), we mention the case of the SQ lattice in the following. The dispersion of the 
SQ lattice is given by E^*^ = —2ti{coskx + cosfcj,) + 4t2COskxCOsky. We take ti=1.0 here. By using 
the above dispersion relations, we obtain the bare Green's function, Go{k, e„) = where 
e„ = 7rr(2n + 1) is the Matsubara frequency and /x is the chemical potential. Ejt is the dispersion 
relation corresponding to each lattice structure. The particle number density n per spin is given 
by " = j^^k,nGo{k,en). 

An effective pairing interaction Vtopt is given by the perturbation expansion up to the third- 
order term with respect to the on-site Coulomb interaction^^"^^) . The origin of the superconduc- 
tivity is investigated by all effective pairing interaction of Vtopt- However, we dare to divide it 
into two parts for convenience, to analysis the role of Vtopt in detail; 

VTOPT(g, k) = FRPA(g, k) + Vvertex(g, k). (2.3) 

The RPA-like term Frpa respects the term included by the Random Phase Approximation (RPA) 
and Vvertex is the vertex correction. The RPA-like term reflects the nature of the simple spin 
fluctuations. The vertex correction term originates from the electron correlation other than the spin 
fluctuations. For the singlet pairing, the RPA-part and the vertex correction part are respectively 
given by 

ySt\Q^ k)=U + U\oiq -k) + 2U\l{q - k), (2.4) 

V^Siq, h) = ^[U'iY: Go{q + k'-k)x (xo(g + k') - Mq + k'))Go{k')] 

k' 

+ U^[J2 Goi-q + k'-k)x ixoi-q + k') - ^o(-g + fc'))Go(A;')]], (2.5) 

k' 

where k indicates k = {k, a;„). The bare susceptibility Xo(<z) and (poiq) are defined respectively by 

Mq) = E Go{k)Go{q + k), (2.6) 

k 

Mq) = ~J2Go{k)Go{q-k). (2.7) 

k 

For the singlet pairing, the Coulomb interaction U connects only the electron lines which have the 
opposite spins. The diagrams for the pairing interaction are shown in Fig. 1. The two external 
lines have the opposite spins. 
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For the triplet pairing, the RPA-like term and the vertex correction correspond to the second- 
order and the third-order terms, respectively. 

V^Jt\<l,k) = -U\o{q-k), (2.8) 

^V™fex*(9, k) = ^[C/'E Go{q -k + k')x {xo{q + k') + Mq + k'))G^{k')] 

k' 

+ C/3[E Go{-q -k + k')x ixoi-q + k') + ^-q + k'))Go{k')]]. (2.9) 
k' 

The diagrams giving the effective interaction for the triplet pairing are shown depend on the 
direction of D-vector. When the two external lines have the parallel spins, the typical diagrams 
are illustrated in Fig. 2. Here, we add the proviso of the diagram in Fig. 1 and Fig. 2. At first 
sight, it seems that the diagram of triplet pairing in Fig. 2 is different from that of singlet pairing 
in Fig. 1. However, the diagram in Fig. 2 is identical with that of Fig. 1.; The first-order diagram 
of T^rI^'^* in Fig. 1 vanishes due to the triplet pairing symmetry. The second and third-order 
diagram of I^r^^^* cancel out each other in the triplet pairing case. Consequently, the residual 
singlet diagram in Fig. 1 agrees with the triplet diagram in Fig. 2. 

An anomalous self-energy is written by the effective interaction V{q, k) and an anomalous Green 
function F{k), as S^(q') = —^TikF{k)V{q,k). At the transition temperature Tc, the value of the 
anomalous self-energy is small and we linearize the Eliashberg equation with respect to F and 
S^, as F{k)^ = \Go{k)\^T,A{k)^ . From these formulae, we obtain the following equation for the 
anomalous self-energy; 

A4(g) = -^J2nq,k)\Goik)f^\ik). (2.10) 

k 

This equation is the linearized Eliashberg equation, which is an eigenequation with an eigenvalue 
A and an eigenvector S^. V{q,k) is given by (2.3), (2.4), (2,5), (2.8) and (2.9). We solve the 
linearized Eliashberg equation on the assumption that has the pairing symmetry shown in 
Table I. We take up all p-, d- and /-wave pairing symmetries. The g-wave is high energy, and the 
possibility of the g-wave pairing is not stronger than that of p-, d- and /-wave pairing in the actual 
superconductivity. Therefore, we investigate the g-wave pairing symmetry in detail. Wc calculate 
the g^yi^x^ -y2y wave pairing as the typical example of the g-wavc pairing symmetry, to obtain the 
full superconducting phase diagrams. The pairing symmetries arc degenerate due to the space 
symmetry in the cubic lattice structure. We show the degree of degeneracy in the cubic lattice 
structure in Table. I. Here, we mention the possibility that the degenerate pairing symmetries are 
mixed below Tc. 

The most dominant pairing symmetry has the largest vahic of the eigenvalues among different 
symmetries. When the eigenvalue A reaches unity, the superconducting state is realized. We 
solve the equation to obtain the dominant state and determine the superconducting transition 
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Table I. Pairing symmetry 



pairing symmetry 


function of symmetry 


degree of degeneracy 




cos(A;2:) — cos(/Cy) 


2 


dxy 


sin(A;a;)sin(A;j^) 


3 


Px 




3 


fxyz 


sin {kx)sm{ky) sin ( kz ) 


1 


fx(3z^—r'^) 


sin(A;^)(2cos(A;2) — cos(A;a;) — cos(A;2^)) 


3 


fz(x'^—y'^) 


sin(A;2)(cos(A;a;) — cos{ky)) 


3 


9xy{x'^—y'^) 


sm{kx)sm{ky) {cos{kx) — cos(/cy)) 
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temperature Tc. The transition temperature Tc is higher than the actual transition temperature, 
because the effect of the normal self-energy is not considered in this calculation. 

§3. Calculated Results 

We divide the first-Brillouin-zonc into 64x64x64 momentum meshes and take Nf = 512 for 
Matsubara frequency. The bandwidth is a necessary region of Matubara frequency a;„ for 
rehable calculations over the low density. The region is covered with the condition; < irTNf. 
To satisfy the condition, we calculate in the region with T >0.004 for the SC lattice structure, 
T >0.006 for the BCC lattice structure and T >0.008 for the FCC lattice structure. In addition to 
this, the condition for U is U <W. 

We obtain the superconducting phase diagram for the electron density n and the next-nearest 
hopping integral ^2- Moreover, we clarify the mechanism of the superconductivity from the eigen- 
value of Eliashberg equation. It depends on the Fermi surface, the bare susceptibility xo and the 
density of state p. Finally, we show the behavior of Tc for various values of U. 

3. 1 Overview for Superconducting Phase Diagram 

In Fig. 3, we show the superconducting phase diagrams determined by the TOPT in the plane 
of the hopping integral t2 and the density n. We figure the most dominant pairing state which has 
the largest value of eigenvalues A. The lattice structures are 2-dimensional SQ and 3-dimensional 
SC, BCC and FCC lattice structures. We can find the following important factors in common with 
the various lattice structures for the SQ, SC, BCC and FCC lattice structures. 

Near the half-filling (n = 0.5), the singlet state becomes dominant by the RPA-like term reflecting 
the simple spin fluctuation process. In the SC lattice, the singlet pairing also exists for t2=-0.5 - 
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0.2 and the intermediate density. 

For the triplet state, the vertex correction plays the vital role for the triplet pairing in the 
intermediate density (0.075 < n < 0.3). On the other hand, in the very low density (n <0.075), 
the triplet pairing state is dominant owing to the simple ferromagnetic spin fluctuations, which 
follows the paramagnon theory^^'^^) The triplet superconductivity can realize in the wide region of 
the density by adding the vertex correction, in contrast to the results^"^) on the basis of only the 
ferromagnetic fluctuations. All phase diagrams in Fig. 3 suggest that the intermediate density is 
necessary for the triplet pairing owing to the vertex correction, in common with the 2-dimensional 
SQ, 3-dimensional SC, BCC and FCC lattice structures. 

Both triplet and singlet superconducting states originate from the general effective interaction 
including both the RPA-like term and the vertex correction. Therefore, the general effective inter- 
action obtained by our perturbation theory gives the unified understanding on the mechanism of 
various superconducting states. In the following subsections, we detail the results for the pairing 
symmetry in the SC, BCC and FCC lattice structures. 

3.2 Simple Cubic Lattice 

^ Superconducting phase diagram 

Fig. 3(b) is the superconducting phase diagram of the SC lattice structure. Near the half-filling, 
the d^2_y2-wave pairing state is dominant. Moreover, the d^2_y2- and dxy-wave pairings also exist 
for t2=-0.5 - 0.2 and the large t2, respectively, from intermediate to high density . 

The p and fxyz-^^^^ pairing states are realized in the intermediate density (0.075 < n < 0.3). 
The triplet state is promoted by the wave number dependence of the vertex correction. In particu- 
lar, the /a;j/2-wave pairing is dominant for the large ^2- Around the /3;y2-pairing region, the p-wave 
pairing is dominant owing to the RPA-like term reflecting the ferromagnetic spin fluctuation. 

<0> Behaviors of eigenvalue and bare susceptibility 

In Fig. 4, we show the n- and t2-dependence of the eigenvalues A obtained by TOPT, the vertex 
correction term and the RPA-like term, respectively, for the p- and d^2_y2-wave pairing states. In 
Fig. 5, we show the bare susceptibility xoiQ^'^n = 0) and the Fermi surface for the various density 
n. From the behaviors of A and xo, we clarify each effect of the RPA-like term VaPA and the vertex 
correction Vvcrtcx- 

The behaviors of A in Fig. 4(a) shows that the d.j.2_y2-wave pairing is induced by the RPA- 
like term near the half-filling. The vertex correction suppresses the d^2_y2-wave pairing. For the 
dx2_y2-wave pairing, xo in Fig. 5(a) shows the spin fluctuations corresponding to a peak of xo at the 
momentum q = (7r,7r,7r). The antifcrromagnctic spin fluctuation originates from a nesting feature 
of the Fermi surface and is reflected by the RPA-like term. Thus, the RPA-like term promotes the 
attractive force for the dj.2_y2--waYe pairing state. 
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For the dxy-wsve pairing, the spin fluctuations at g = (0, it, 0), (tt, 0, 0), (0, 0, tt) are reflected by 
the RPA-hke term in Fig. 5(b). Thus, the dxy-wave pairing is dominant for the large next-nearest 
hopping integral t2, which gives a strong nesting feature. The vertex correction also suppresses the 
dxy-wscve pairing. 

In Fig. 4(b), 4(c) and 4(d), we show the behaviors of A for the p- and fxyz-^^^^ pairing states. 
When the solution of the Eliashberg equation does not converge, we set A zero. For the intermediate 
density, the p- or fxyz-^^'^^ pairing state has the largest eigenvalue owing to the vertex correction. 

The superconducting mechanism due to the vertex correction is explained in the following. The 
vertex correction has a unique wave number dependence. The wave number dependence of the 
vertex correction is not owing to the simple spin fluctuation process, because any character of the 
spin fluctuation does not clearly appear in Fig. 5(c) and 5(d). The characteristic wave number 
dependence of the vertex correction is induced by the general particle-particle iteration due to the 
Coulomb interaction. The wave number dependence induces the attractive force between the elec- 
trons near the Fermi surface. The RPA-like term suppresses the eigenvalue for the p- and f-wsye 
pairing states. 

<0> Behaviors of eigenvalue and density of states 

In Fig. 4(b), A for the p-wave pairing is large from n=0.175 to 0.275. This is because the density 
of states p is large at Fermi level from n=0.175 to 0.275 in Fig. 6(a). In Fig. 4(c), A for the 
p-wave pairing is small as t2 is small. The decrease originates from the result that p at Fermi level 
decreases as t2 becomes small, which is shown in Fig. 6(b). 

3.3 Body- centered Cubic Lattice 

<0> Superconducting phase diagram 

Fig. 3(c) is the superconducting phase diagram of the BCC lattice structure. In the high density 
near the half-filling, the dxy-wa,ve pairing is induced by the RPA-like term which reflects the simple 
spin fluctuations. In the intermediate density (0.075 < n < 0.3), the p-wave pairing is promoted 
by the wave number dependence of the vertex correction. The /- and ^-wave pairing states do not 
appear in this case. 

<0> Behaviors of eigenvalue and bare susceptibility 

In Fig. 7(a) and 7(b), we show the n- and t2-dcpcndcncc of the eigenvalue A obtained for 
the p-wave pairing state. The behavior of A is similar to the case of the SC lattice structure. 
The p-wave pairing state is dominant owing to the vertex correction in the intermediate density 
(0.075 < n < 0.3). The RPA-like term suppresses the p-wavc pairing. On the other hand, the 
dxy-wavc is dominant owing to the RPA-like term near the half-filling, while the vertex correction 
suppresses the dxy-wave pairing. This mechanism for the p- and d-wave pairing states are similar 
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to that in the case of the SC lattice. 

In Fig. 8, we show the bare susceptibihty Xo{Q,i^n = 0) and the Fermi surfaces for n=0.45 
and 0.15 at t2=0.1, respectively. In Fig. 8(a), the antiferromagnetic spin fluctuations at 
q = (tt, 7r,7r), (tt, 0, 0) and (0,7r,0) encourage the dxy-wave pairing state. The antiferromagnetic 
spin fluctuations are induced by the nesting feature of the Fermi surface, such as the dxy-wsve 
pairing case of the SC lattice. In Fig. 8(b), Xo{Q,(^n = 0) does not have any character of the spin 
fluctuation but the vertex correction induces the p-wave pairing. 

Behaviors of eigenvalue and density of states 

We show the n- and t2-dependence of the density of states p in Fig. 9. A is large from n=0.4 to 
0.475 in Fig. 7(a), because the density of states p in Fig. 9(a) is large at the Fermi level. The A 
in Fig. 7(b) becomes small as t2 is small. This is because that the density of states at Fermi level 
decreases as t2 becomes small, which is shown in Fig. 9(b). 

3.4 Face- centered Cubic Lattice 

<^ Superconducting phase diagram 

In Fig. 3(d), we show the superconducting phase diagram for the FCC lattice structure. Near 
the half-filling, the dxy- or g-wave pairing state is dominant by the RPA-like term. The p-wave 
pairing state is promoted by the vertex correction in the intermediate density (0.075 < n < 0.3). 
In contrast with this, the p-wave pairing is induced by the RPA-like term for the low hole density 
(n > 0.7). 

<0> Behaviors of eigenvalue and bare susceptibility 

In Fig. 10, we show the n- and t2-dependcncc of the eigenvalue A for the p-w&ve pairing state. 
In the intermediate density (0.075 < n < 0.3), the p-wave pairing state is dominant. The RPA-like 
term suppresses the p-wave pairing. On the other hand, the RPA-like term encourages the triplet 
pairing in the low hole or low electron density (n > 0.7 or n < 0.075). Thus, the ferromagnetic fluc- 
tuations play the important role in the triplet superconductivity. The vertex correction suppresses 
the p-wave pairing. 

In Fig. 11, the bare susceptibility xo{Qi'-^n = 0) and the Fermi surface arc shown. Near the 
half-filling, there are the spin fluctuations corresponding to a peak of xo stt a certain momentum 
vector in Fig. 11(a), 11(b) and 11(c). The spin fluctuations arc reflected by the RPA-like term and 
the term gives the advantage for the d^2_y2-, dxy- and g-wavc pairing states. The vertex correction 
suppress the singlet pairing states. The spin fluctuations are induced by a feature of the Fermi 
surface. The spin fluctuations at g = (7r/2,7r/2, k^) encourage the d^2_y2-wavc pairing state in Fig. 
11(a). On the other hand, the spin fluctuations at q = (0, 0, vr), (0, 7r,0), (tt, 0, 0) encourage the 
dxy-wave pairing state in Fig. 11(b). 
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For the intermediate density (0.075 < n < 0.3) in Fig. 11(d), the feature of a certain spin 
fluctuation does not clearly appear. The wave number dependence of the vertex correction induce 
the p-wave pairing. 

For the low hole or low electron density, xoiQ^^n = 0) has the ferromagnetic spin fluctuation, 
which exhibits a peak of xo near q = (0,0,0) in Fig. 11(e). The ferromagnetic spin fluctuation 
gives the advantage of the p-wave pairing state. The situation of the low hole density is connected 
with the characteristic density of states in the FCC lattice, which also induces the ferromagnetism. 
The detail of explanation is the following. 

Behaviors of eigenvalue and density of states 

The density of states p is large at the Fermi level as shown in Fig. 12(a) in the low hole density 
(n = 0.84). This is the particular nature of the density of states in the FCC lattice structure. On 
the basis of the particular band structure, the ferromagnetism of Ni is explained by Kanamori^^. 
In the low hole density, the p-wave pairing is dominant owing to the ferromagnetic fluctuations, 
which agrees with the Arita's results^^. 

However, the triplet pairing induced by the ferromagnetic spin fluctuation does not have the 
large eigenvalue in the Eliashberg equation. This fact indicates that it is not easy for the simple 
ferromagnetic spin fluctuation to introduce the triplet pairing state. Moreover, the eigenvalue is 
dominant due to the ferromagnetic spin fluctuation and is suppressed by the vertex correction, as 
the electron correlation becomes strong. Therefore, the superconducting transition temperature is 
very low for the realistic superconductor to exist owing to the ferromagnetic spin fluctuation. 

In Fig. 10(c), A also decreases as t2 becomes small. The decrease originates from the fact that p 
at Fermi level decreases as t2 becomes small as shown in Fig. 12(b). 

3.5 Dependence of Eigenvalue on U and effect of vertex correction 

In this subsection, we certify that the triplet pairing induced by the vertex correction exists in 
the wide region of U. 

We show the f7-dependence of the eigenvalue A obtained on the basis of TOPT in Fig. 13. We 
focus the density n at the crosspoint between the p-wave and the d-wave pairing. The p-wave 
pairing is dominated by the vertex correction. The crosspoint shifts to small n as U becomes small 
from 7.0 to 2.0. The bandwidth W in the SC lattice structure equals 12. Thus, the density region 
of the p-wave pairing state becomes narrow, while the density region of the d-wave spreads. The 
density n at the crosspoint changes from 0.22 to 0.32 in Fig. 13. The change is not so large, and 
the phase diagram in Fig. 3(b) does not change so much for the value of U, as seen in our previous 
resultB of the 2-dimensional SQ lattice. Therefore, the p-wave due to the vertex correction exists 
in the intermediate density in the wide region of U. For the small U /W such as U = 2.0 in Fig. 
13(a), the effect of the higher-order interactions is weaker than that of the second- and third-orders. 
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(For U < 2.0, we do not calculate A, because A is too small to obtain a reliable calculation.) Even 
if we take into account of the effect of the higher-order interactions, the p-wave induced by the 
vertex correction seems to be dominant in the intermediate density. This fact gives reliability to 
the result of TOPT. 

3.6 Transition Temperature T^ and dimensionality 

We discuss the effect of the dimensionality on Tc- In Fig. 14, we show the n-dependence of the 
superconducting transition temperature T^ calculated by TOPT for the p- and d^,2_j^2-wave pairing 
states. The parameters are t2=0.1 and C/=10.0. The unit of energy is the hopping transfer ti. T^ in 
the 3-dimensional SC lattice is one-order lower than Tc in the 2-dimensional SQ lattice structure.0) 
Therefore, both singlet and triplet pairing states are in favor of the 2-dimensional system such as 
the layered system rather than the 3-dimensional system. 

§4. Summary and Conclusion 

We have studied the superconductivity in the Hubbard model by the third-order perturbation 
theory. We conclude by pointing out the main factors in common with the 2- and 3-dimensional 
various lattice structures; the SQ, SC, BCC and FCC lattice structures. 

(I) Study on the basis of effective interaction including vertex correction 
Singlet pairing 

For the singlet pairing, the spin fluctuation is the important factor. Near the half-filling or the 
intermediate density, the singlet superconductivity is realized by the spin fluctuation, such as the 
antiferromagnetic spin fluctuation. The vertex correction suppresses the singlet superconductivity 
induced by the spin fluctuation. 

In the 3-dimensional system, the effect of the magnetic order is strong near the half-filling. There- 
fore, it seems to be difficult to realize the unconventional singlet superconductivity. On the other 
hand, in the heavy fermion system, the magnetic instability is usually suppressed by the Kondo 
effect.lli' In this case, the singlet pairing is dominant around the half-filling case in the supercon- 
ducting phase diagram of the SC lattice. The cubic heavy fermion superconductor Celns may be 
a good example. 

Triplet pairing 

By taking account of the vertex correction, the triplet pairing can be realized in the wide region 
of the electron density. The triplet pairing is induced by the vertex correction in the intermediate 
density. 

On the other hand, the p-wave pairing is also induced by the ferromagnetic spin fluctuation 
process for the low hole or low electron density. In particular, the case of the low hole density in 
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the FCC lattice has a similar electronic structure to that of the ferromagnetism of Ni and nearly 
ferromagnetic Pd. The density of states in the FCC lattice has the narrow and high density peak 
at the Fermi level for the low hole density. The narrow and high density of states gives rise to 
the ferromagnetic spin fluctuation, which induces the triplet pairing state. This case of FCC lat- 
tice is the particular case in contrast with the SQ, SC and BCC lattice structures. The triplet 
pairing induced by the ferromagnetic spin fluctuation does not have the large eigenvalue in the 
Eliashberg equation. Moreover, our theory shows that the eigenvalue owing to the ferromagnetic 
spin fluctuation is suppressed by the vertex correction as the electron correlation becomes strong. 
Therefore, the superconducting transition temperature is too low to realize the superconductivity 
in the strong correlation system. This fact indicates that it is not easy for the simple ferromagnetic 
spin fluctuation to introduce the triplet pairing state. The result explain the fact that the super- 
conductiviting state has not been actually found in the Ni and Pd metals, although there exist the 
strong ferromagnetic spin fluctuation. 

On the other hand, in the very low electron density, the triplet state is also dominant owing to 
the ferromagnetic spin fluctuation. However, the very low density does not apply to the electron 
density of the realistic material. In the strong coupling system, the realistic materials usually have 
the intermediate or high density. The superconducting phase diagram (Fig. 3) indicates that the 
triplet pairing induced by the simple ferromagnetic spin fluctuation is confined to the smaller region 
than that induced by the vertex correction in the intermediate density. Therefore, it seems that the 
triplet pairing is difficult to be induced by the simple ferromagnetic spin fluctuation theory such as 
the paramagnon theory. We think it is natural that the triplet pairing is induced by the effective 
interaction including the vertex correction. 

^ All effective interaction including vertex correction 

The effect of the vertex correction is generally opposite to the effect of the RPA-like term. When 
the vertex correction mainly gives the advantage to the triplet pairing, the RPA-like term suppresses 
the triplet pairing. When the RPA-like term reflecting the antiferromagnetic (ferromagnetic) spin 
fluctuations induces the singlet (triplet) pairing state, the vertex correction suppresses the pairing. 

When we investigate the unconventional superconductivity induced by the strong Coulomb cor- 
relation, we should consider the effective interaction as the origin of the unconventional supercon- 
ductivity. In this case, it seems necessary to study not only the spin fluctuations but also the wave 
number dependence of the vertex correction. By analyzing the general wave number dependence of 
the quasi particle interaction on the basis of Fermi liquid theory, it is possible to derive generally 
both the singlet and triplet superconductivity. 

(II) Dimensionality and actual superconductors 
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Tc in the 3-dimensional system is one-order lower than that in the 2-dimensional SQ lattice. 
This result indicates that the superconductivity is in favor of the 2-dimensional system such as 
the layered structure. The many unconventional superconductors exist in the 3-dimensional heavy 
fermion system. The unconventional superconductivities in the 3-dimensional system seems to exist 
in Celn30) and the filled skutterudite PrOs4Sbi2i3) . 

Celns contains the SC lattice constructed by Ce atoms. Celna is the unconventional supercon- 
ductoiH* and the 3-dimensional system composed of Ce 4/ electrons. © The topology of the main 



3-dimensional Fermi surfacaHs' is similar to the d^2 _y2-wa,ve pairing case of Fig. 5(a) in the SC 
lattice. Near the half-filling in the 3-dimensional SC lattice, the antiferromagnetic order is strong 
near the half-filling in the 3-dimensional system. In the heavy fermion systems, the antiferromag- 
netic order is not so strong, and the superconductivity is also possible in the 3-dimensional SC 
lattice. Therefore, the superconductivity of Celna may be explained by the d-wave pairing state in 
SC lattice. Our result of the SC lattice has suggested that the d^2_y2-wa.ve pairing is induced, in 
the case that the antiferromagnetic spin fluctuation exits near Q=(7r/2, 7r/2, 7r/2) at Ce site. The 
low transition temperature Tc ~ 0.2 K of Celna matches Tc of the 3-dimensional system lower than 
that of the 2-dimensional system, which is shown in our result. Actually estimating the real Tc in 
Celn3, we consider the effect of the Fermi energy which is reduced by the renomarization factor 
z ^ 0.05. 

The filled skutterudite PrOs4Sbi2 has the 3-dimensional Fermi surfacJH^ of the cubic lattice. In 
the study of the Sb-NQR of the PrOs4Sbi2,ii* the dependence of 1/Ti on T does not have a clear 
coherence peak, while it has the exponential behavior in the low temperature. The fact may indicate 
the possibility of the unconventional superconductivity with the full superconducting gap such as 
the BW state of ^He-superfluidity .S^ The superconducting gap symmetry might change such as 
UPt3,0'* because the specific heat has the behavior of the double transition .0^ This superconductor 
may be explained as the triplet pairing state owing to the effective interaction including the vertex 
correction. According to the present theory shown in Fig. 3, odd parity or even parity states 
appear mainly depending the electron density n and the feature of Fermi surface. 

Applying this result to the 3-dimensional Fermion system, we should reconsider the origin of the 
superfluidity in ^He. To make clear the origin of the superfluidity in ^He@ as the 3-dimensional 
system, it seems necessary to study not only the simple spin fluctuations but also the effective 
interaction including the wave number dependence of the vertex correction. This study is partly 
given in the papers^^'^'^. 

(Ill) Summary 

Arita et al. investigated the pairing instability and the magnetic behavior in the 2- and 3- 
dimensional cubic lattice system on the basis of the simple spin fluctuations by the FLEX approx- 
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imation.EP For the superconductivity owing to the antiferromagnetic (ferromagnetic) spin fluctu- 
ation, our result agrees with their's results. Therefore, TOPT is the reliable theory covering the 
theory of the spin fluctuation. 

The triplet pairing state is induced by the effective interaction including the vertex correction 
rather than that by the simple ferromagnetic fluctuation. The effective interaction including the 
vertex correction plays the important role for the triplet superconductivity in the 3-dimensional 
system as well as the 2-dimensional system-'^-^'-'^^). When we investigate the unconventional super- 
conductivity in the strong correlated system, the mechanism should be considered owing to not 
only the simple spin fluctuation but also the effective interaction including the vertex correction. 
A part of this study is given in the papers.ii^ 
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Fig. 1. Diagrams for the effective interaction of singlet pairing within the third-order perturbation with respect to 
U. The solid line is the bare Green's function Go- The broken line is the Coulomb interaction U. The broken line of 

U connects only solid lines possessing opposite spins. The two external lines have the opposite spins. The effective 
interaction is divided into the RPA-like part and the vertex correction. The latter begins with the third-order 
terms. 



Fig. 2. Diagrams for the effective interaction of the triplet pairing within the third-order perturbation. The two 
external lines have the parallel spins. The RPA-like part and the vertex correction of the pairing interaction are 
given by only the second-order and the third-order terms, respectively. 
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Fig. 3. The superconducting phase diagrams for the dominant pairing symmetry. t2 and n are the next-nearest 
hopping integral and the density, respectively. The half-filling density corresponds to n=0.5. 
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Fig. 4. The n- and t2-dependence of the eigenvalue A obtained by TOPT, vertex and RPA-like terms for the p-, 
d^2_y2- and /^j/z-wave pairing states in the SC lattice structure. 
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Fig. 5. The bare susceptibility xo(5, = 0) and Fermi surface in a quarter first-Brillouin-zone of the SC lattice 
structure, (a) and (b) are in the case that the d-wave pairing states are induced by the spin fluctuation, (c) and 
(d) are in the case that p- and /-wave pairing states are dominant owing to the Vertex correction. 
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Fig. 6. The n- and i2-dependence of the density of states at T—0 and U=0 for the SC lattice structure. The 
parameters are fixed t2=0.1 for (a) and n=0.225 for (b), respectively. 



19 



Fig. 7. (a) The n-dependence of the eigenvalue A for the p- and d^2_y2-vfaye pairing states for the BCC lattice 
structure. The parameters of (a) are fixed at t2=0.1, U=6.0 and T=0.006. (b) The t2-dependence of the eigenvalue 
A for the p- and d^2_j,2-wave pairing states for the BCC lattice structure. The parameters are fixed at n=0.15 and 
T=0.008. 
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Fig. 8. The bare susceptibility xo(?,w„ = 0) and Fermi surface for the BCC lattice structure for the BCC lattice 
structure, (a) is in the case that the dxy-pairing state is induced by the spin fluctuation, (b) is in the case that the 
p-wave pairing is encouraged by the Vertex correction. 



21 



Fig. 9. The n- and t2-dcpcndcncc of the density state at T=0 for the BCC lattice structure. The parameters are 
fixed n=0.175 and t2=0-l, respectively. 
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Fig. 10. (a) The n-dependence of the eigenvalue A for the p-wave pairing state for the FCC lattice structure. The 

parameter is fixed at t2=0.1. (b), (c) The f2-dependence of the eigenvalue A for the p-wave pairing state for the 
FCC lattice structure. The parameters are fixed at n=0.15 and 0.84, respectively. 
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Fig. 11. The bare susceptibility xo{'l,'^n ~ 0) and Fermi surface for the FCC lattice structure, (a), (b) and (c) 
exhibit xoil, = 0) in the cases of the d- and g-wave pairing states due to the RPA-like term, (d) corresponds to 
the p-wave pairing which is encouraged by the Vertex correction term, (e) corresponds to the p-wave pairing owing 
to the RPA-like term. 
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Fig. 12. The n- and t2-dcpcndence of the density of states at T=0 for the FCC lattice structure. The parameters 
are fixed ^2=0.1 and n=0.15, respectively. 
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Fig. 13. The f/-dependence of the eigenvalue A obtained by TOPT, vertex and RPA-hke terms for the p-, d^2_y2- 
wave pairing states in the SC lattice structure. 
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Fig. 14. The n-dependence of the superconducting transition temperature Tc for the p- and d^2_y2-wauve pairing 
states for the SC lattice structure. Tc is calculated by TOPT. The parameters are t2=0.1, U=W.O. The unit of 
energy is the hopping transfer ti. 
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